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Abstract 

In realistic unified models involving so-called 5O(10)-inspired patterns of Dirac 
and heavy right-handed (RH) neutrino masses, the lightest right-handed neu¬ 
trino N\ is too light to yield successful thermal leptogenesis, barring highly fine 
tuned solutions, while the second heaviest right-handed neutrino N 2 is typically in 
the correct mass range. We show that flavour coupling effects in the Boltzmann 
equations may be crucial to the success of such IV 2 dominated leptogenesis, by 
helping to ensure that the flavour asymmetries produced at the IV 2 scale survive 
N\ washout. To illustrate these effects we focus on IV 2 dominated leptogenesis 
in an existing model, the A to Z of flavour with Pati-Salam, where the neutrino 
Dirac mass matrix may be equal to an up-type quark mass matrix and has a 
particular constrained structure. The numerical results, supported by analytical 
insight, show that in order to achieve successful N 2 leptogenesis, consistent with 
neutrino phenomenology, requires a “flavour swap scenario” together with a less 
hierarchical pattern of RH neutrino masses than naively expected, at the expense 
of some mild fine-tuning. In the considered A to Z model neutrino masses are 
predicted to be normal ordered, with an atmospheric neutrino mixing angle well 
into the second octant and the Dirac phase 5 — 20°, a set of predictions that 
will be tested in the next years in neutrino oscillation experiments. Flavour cou¬ 
pling effects may be relevant for other S'O(10)-inspired unified models where N 2 
leptogenesis is necessary. 
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1 Introduction 


It would be certainly desirable to extend the SM, realising a unified picture able to solve 
the flavour problem, explaining masses and mixing parameters of quarks and leptons, 
and at the same time to provide solution to the cosmological puzzles. In particular 
leptogenesis [1] is an attractive way to explain the matter-antimatter asymmetry of the 
Universe, that in terms of the baryon-to-photon number ratio is given by [2] 

'^MB = ( 6 .1 ± 0 .1) X 10 "1° . (1) 

Leptogenesis is a cosmological application of the see-saw mechanism, a minimal exten¬ 
sion of the SM able to explain neutrino masses and mixing [3]. Despite an intense 
activity on various aspects of leptogenesis, there are not many definite realistic unified 
models that have been shown to lead to successful leptogenesis while explaining fermion 
masses, mixing and CP violation, although in the case of SU (5) it is certainly possible 
to achieve successful N\ leptogenesis (for a recent example see e.g. [4] which uses the 
sequential dominance results for N\ leptogenesis discussed in [5]). 

In this paper we are interested in N 2 leptogenesis in so-called S'O(10)-inspired models 
with type I seesaw. For definiteness, we investigate the possibility that the A to Z 
Pati-Salam model proposed in [6] (see also [7]) can not only describe neutrino masses 
and mixing but also attain the correct value of the matter-antimatter asymmetry with 
leptogenesis. The right-handed (RH) neutrino mass spectrum in this model is very 
hierarchical, typical of S'O(10)-inspired models (assuming type I seesaw). In this way 
the lightest RH neutrino Ni is too light to generate a sizeable asymmetry [8] while on 
the other hand the next-to-lightest RH neutrino is heavy enough to be able potentially 
to generate the correct asymmetry realising the so called N 2 dominated scenario [9, 10, 
11, 12] or simply N 2 leptogenesis. However in general it is non-trivial to be able to find 
a set of values of the parameters of the model for the lightest RH neutrino wash-out 
to be negligible, in such a way that the asymmetry generated by the next-to-lightest 
neutrinos survives, while at the same time producing values of the neutrino parameters 
compatible with the experimental results. 

The non-trivial requirements of N 2 leptogenesis may be somewhat relaxed by taking 
into account additional flavour coupling effects in the Boltzmann equations, or more 
generally the kinetic mixing terms, that may transmit part of the initially produced 
flavour asymmetry from one particular flavour to other flavours. This increases the 
chances that the asymmetry generated in a particular flavour at the N 2 scale may 
survive the Afi washout. The additional flavour coupling effects, which are not usually 
considered in the literature, arise mainly from the fact that lepton asymmetry produced 
at the N 2 scale into left-handed lepton doublets is also accompanied by hyper-charge 
asymmetry into Higgs asymmetry, giving the dominant effect, baryons asymmetries 
into quarks and lepton asymmetries into right-handed charged leptons, which are then 
transmitted to the other flavours which couple to these particles [13]. This results in 
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new contributions to the total asymmetry that in a traditional N ]-dominated scenario 
would give just some small corrections (0(10%)) [14]. However, in a ^-dominated 
scenario they can become dominant if, contrarily to the usual dominant terms produced 
in a flavour that is strongly washed by the lightest RH neutrinos, they are produced 
into a flavour that escapes the lightest RH neutrino wash-out. In this way the so called 
“flavour swap scenario” is realised [15]. 

In the A to Z of Flavour with Pati-Salam [6], all of these features are exemplified. In 
particular, the Ni scale is too light to generate asymmetries thermally, while the flavour 
asymmetry produced at the N 2 scale, namely the r flavour, is effectively washed out at 
the N\ scale for all ranges of parameters consistent with the experimentally acceptable 
neutrino masses and mixings. In this example, flavour coupling effects come to the 
rescue, effectively transmitting part of the r asymmetry into (a linear combination 
of) electron and muon asymmetries at the N 2 scale, completely (electron) or partly 
(muon) surviving washout at the Ni scale. These features necessarily arise from the 
rigid structure of the Yukawa and Majorana matrices enforced by the model, leading to 
fairly precise predictions for PMNS parameters which can be compared to experiment. 
It has been already shown that the model can simultaneously fit both lepton and quark 
parameters. Here for simplicity we will focus on the lcptonic sector and consider non- 
supersymmetric leptogenesis. 

The layout of the remainder of the paper is as follows. In Section 2 we review neutrino 
masses in the A to Z model. In Section 3 we fit the neutrino parameters without imposing 
successful leptogenesis. In Section 4 we show how within a traditional calculation of 
the asymmetry ignoring flavour coupling effects one cannot find any good fit both to 
neutrino parameters and to the measured value of the asymmetry simultaneously. In 
Section 5 we show how flavour coupling rescues the model. In Section 6 we provide an 
analytical insight on the found solution within more general context of SO (10)-inspired 
leptogenesis highlighting different aspects both of the specific solution within the A 
to Z model and more generally of S'O(10)-inspired models. In Section 7 we draw our 
conclusions. 


2 Neutrino masses in the A to Z model 

The lowest order lepton Yukawa matrices (in LR convention) and heavy Majorana mass 
matrix Mr are predicted by the model just below the high energy Pati-Salam breaking 
scale ~ few x 10 16 GeV. The charged lepton Yukawa matrix of the model is diagonal to 
excellent approximation, and the neutrino Lagrangian in this basis is given by, 

-C = N l Y'"N' r + Nj M' R N' R + H.c. (2) 
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with the neutrino Yukawa and Majorana matrices [6], 

/ 0 b e ~ i3 ”/ 5 0 \ /M' u e 2 * 0 M’ 13 e*\ 

Y' u = ae -* 37r / 5 4fc-i37r/5 0 ] , M' r — I 0 Ml n e* 0 , (3) 

\ae~ i3lT/5 2 be~ i3n/5 ce*) \ M[ 3 A 5 6 M’ 33 ) 

where M' u , M[ 3 , M 22 , M 33 are positive and real. Note that for simplicity it was also 
assumed in [6] that M[ 3 = 0 and 0 = 0, while here we shall allow these parameters 
(generally present in the model) to take non-zero values. Generally the model predicts 
M[ x < M 22 M 33 and M[ 3 ~ M 22 . Without any tuning of parameters, the model 

predicts the typical values, 

M’ u ~ 10 5 GeV, M 22 ~ M[ 3 ~ 10 10 GeV, M' 3 ~ 10 15 GeV, (4) 

while 

a : b : c ~ m up : m charm : m top ~ 10~ 6 : 10” 3 : 1. (5) 

However, successful leptogenesis requires some fine-tuning, leading to a more compressed 
spectrum of right-handed neutrino masses, as we shall see. Also £ is chosen to be one 
of the complex fifth roots of unity: £ = 0, ±27r/5, ±47 t/ 5. We shall consider all cases. 

We perform a unitary transformation Ur to the flavour basis where the right-handed 
Majorana mass matrix is diagonal with real, positive eigenvalues, 

£ = N L Y’‘'U R U' R N' R + NjU' R UlM' R U R U) t N' R + H.c. (6) 

= N l Y-N r + N t r M r N r + H.c. (7) 

where 

Nr = [/£> N' r , Y" = Y v U R , Mr = Ur M' r U r = diag(M 1 , M 2 , M 3 ). (8) 

We parametrise the unitary matrix, assuming small angle rotations, approximately as 

/ 1 0 i?i 3 e^ 13 \ /e^ 11 0 0 \ 

U R = ( 0 1 0 0 e^ 22 0 , (9) 

0 1 / \ 0 0 e^ 33 / 

where the first matrix factor diagonalises the right-handed neutrino mass matrix and the 
second factor ensures that the right-handed neutrino masses M* are real and positive. 
Thus the form of the neutrino Yukawa matrix in the flavour basis is Y u = Y v Ur, 

/ 0 6e —7.37r/5 e 7^ 22 Q \ 

Y u ~ I ae- i:i7T/5 e l4> " 4&e -;3*-/5 e # 22 a _r 13 e <(03s-*-37r/5) ,( 10 ) 
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The parameters of Ur are determined from the requirement that 


Mr — U t r M' r U r — Dm = diag(Mi, M 2 , M 3 ). 

In particular, the matrix elements of the diagonal Mr must satisfy: 


(ii) 


{M r ) n 

= Mi 

(Mr) 13 

= 0 « 

(Mr ) 22 

= m 2 

(Mr) 33 

= M 3 


( 12 ) 


From these conditions (remembering that M, are real and positive) we find: 


R 


13 


ML 1 

Jjr, 013 = -z, 011 = -C+7^Arg 


leading to: 


Mi 


Mo 


33 


M [| — 2M ( 3 /? I3 + M' i:i R\ 3 ~ 


M' 


M ' 2 


13 


11 AT' 3J 


, 022 « -^2, 033 ~ 0 ,(13) 


11 ML ’ 


*33 


m' 2 , 


M 3 ~ M 33 + M n i ? 2 3 + 2M 13 i?i 3 . 

Thus the neutrino Yukawa matrix in the flavour basis is given by, 
y" % 


0 frg-* (£/2+3 tt/5) q 

a g-*(?+37r/5) 45e _ * (^/ 2+37r / 5 ) a _R 13 e -*(?+37r/5) 

k a e _i (€+3 tt/5)( 1 _ 7 ) 26e- iK/2 +3^/5) c e # 


where 


7 « (c/a) J Ri 3 e i( ^ +37r/5) . 


(14) 


(15) 


(16) 


The Dirac neutrino masses, the eigenvalues of the neutrino Dirac mass matrix m r J = 
v Y rr , are given by 


rn^ x ~ av 


/y/vf, m - ~ 


£2 ~ vTr b 


, D 


(17) 


where the approximations are quite precise given the working assumptions leading to 
the typical values in Eq. (5). 
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For the up-type quark matrix we have two options, 


/ 0 fee"* 3t/5 0 \ 

Y u = ae~ i3n / b Ybe~ i3n ' b 0 CASE A (18) 

y ae -^/5 2 be-* 3 ** 5 ce i<f> J 

and 

/ 0 | be~ i3n / 5 0 \ 

Y u = ae~ i:i7r ^ | b e ~ i3n / 5 0 CASE B. (19) 

\ae~ i3nlb § 6e" i37r/5 3 ce i4> ) 

These are simply related to the neutrino Yukawa matrix in the original basis, Y >v , by 
Clebsch relations for CASE B, while we have simply Y u = Y v for CASE A (in the basis 
of Eq.(3)) which is just the minimal S'O(10)-like expectation that the Dirac neutrino 
mass matrix is identically equal to the up-type quark mass matrix. Notice that in the A 
to Z model there are “texture zeroes” in the (1,1), (1,3) and (2,3) entries of the Yukawa 
matrices above that will play a role in leptogenesis considerations. 

The Dirac neutrino masses are simply related to the up-type quark masses, depending 
on the choice of model, 

TTl^i m up , 1 ^ 1/2 TTTcharm; ^Y /3 tTltop CASE A (20) 

and 

"7 ~ m up , m ” 2 « 3m charm , m® 3 « ^ m top CASE B. (21) 

In this way the (real) parameters a, b and c are simply determined from the values of 
the up quark masses at the grand-unified scale 1 so that the (real) input parameters are: 

M [,, M( 3 , M 22 , M 33 and phase 0. The phase £ is restricted to be one of the complex fifth 
roots of unity: £ = 0, ±27 t/5, ± 47 t/ 5 . 

Seesaw mechanism in the flavour basis 

Using the see-saw formula 

m v = —v 2 Y v M^ 1 Y uT (22) 

in the flavour basis we And the neutrino mass matrix m u , 

/0 0 0 \ (\ 4 2\ (0 0 0\ 

m v « m a e i4,a 0 1 (1 - 7 ) +m b e i<t>b 4 16 8 +m c e 1 ^ 0 0 0 ,(23) 

\o (1-7) (1-7)7 V 2 8 V \° 0 l ) 

1 We use the values m toP = 100 GeV, TO c harm = 400 MeV and ?n up = 1 MeV. [16] 
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(24) 


a 2 v 2 


rrin = 


Mi 


, m b = 


b 2 v 2 
Mft 


m r = 


c 2 v 2 

mT 


= -2(£ + 3tt/ 5), ft = -2(£/2 + 3 tt/ 5), ft = 2ft 


(25) 


From Eqs.(17),(24) and (20) or (21) we find, 


m n 


m; n 

17- up 


Mi ’ 


m h ~ (9) 


m; 


charm 


m 


17M 2 


rn r 


top 


(9 )M 3 


(26) 


where the factors in brackets apply to CASE B, and these factors are simply unity for 
CASE A. Thus the three right-handed neutrino masses Mi, M 2 , M 3 may be determined 
from m a , m b , m c . 

The Majorana neutrino mass matrix m u , defined by C v = — |m"I7 l^£ + h- c -> i s 
diagonalised by 


(m i 0 0 \ 

U UL rn v Ul=\ 0 m 2 0 . (27) 

\ 0 0 m 3 ) 

The PMNS matrix is then given by 

CpMNS = Ue h Ul h . (28) 

We use a standard parameterization Cpmns = R 23 U 13 R 12 P in terms of s tJ = sinftft, 
Cij = cos(%), the Dirac CP violating phase 5 l and further Majorana phases contained 
in P = diag(e*~2", e*^, 1). The standard PDG parameterization [17] differs slightly due 
to the definition of Majorana phases which are by given by Ppdg — diag(l, ). 

Evidently the PDG Majorana phases are related to those in our convention by a 2 i = 
/ft — ft and a 3 i = —ft, after an overall unphysical phase is absorbed by Lft. 

For example, using the input parameters £ = Att/5, 0 = 7 = 0 and 

m a = 0.034 eV, m b = 0.002 eV,m c = 0.002 eV, (29) 

we find for the CASE A the following values for the physical neutrino masses, 

m 1 = 0.00034 eV, m 2 = 0.0086 eV, m 3 = 0.050 eV. (30) 

and the lepton mixing parameters, 

0i2 = 34.0°, 0i3 = 9.1°, ft 3 = 39.7°, 5 = 260°, ft = 321°, ft = 75° . (31) 
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3 Fitting neutrino parameters 


In this Section we perform a quantitative numerical analysis. We randomly generate a 
set of values of the five (continuous) parameters of the model (y = M[ x , M ! n , M 33 , M[ 3 , (ft) 
for each of the five choices of the discrete parameter £ and from these we calculate the 
nine low energy neutrino parameters in m u . A recent global analysis [18] finds the values 
shown below (with la errors). 

The solar mass squared difference: 

A/U 32 = ml — m\ = 7.50±o;;l7 x 10 _ 5 eV 2 . (32) 


The atmospheric mass squared difference, respectively for normal ordering (NO) and for 
inverted ordering ( 10 ): 


Am^ = m 3 — m i = 2 - 457 ^o;q 47 x 10 3 eV 2 and 

A/773 2 = 777 2 — 7772 = X 10~ 3 eV 2 . 

The solar mixing angle: 

n _ qq aoo +0.78° 

(712 — 00.40 _q 750 . 

The reactor mixing angle: 


The atmospheric mixing angle, respectively for NO and for 10: 

9-23 = 42.3°and 9 23 = 49.5° . 

The Dirac phase, respectively for NO and for 10: 

6 = 306V?’° and 6 = 254VS . 


(33) 

(34) 

(35) 

(36) 

(37) 


We see that, for the atmospheric mixing angle, within la, there are currently two so¬ 
lutions, one in the first octant for NO and one in the second octant for 10, though the 
alternative choice of octant is in both cases only disfavoured at only ~ 1.4 a. Note also 
that, within 3cr, the Dirac phase 6 can have any value. 

The four parameters Am]^, At 77 | 3 , 9 12 and 9\ 3 , are measured quite accurately and 
precisely and their distributions are very well approximated by Gaussian distributions. 
On the other hand the atmospheric mixing angle is not only much less precisely measured 
and not Gaussianly distributed, but is also affected by larger systematic uncertainties 
and in particular different global analyses exclude maximal mixing with different statis¬ 
tical significance, in any case below 2 a, so that its determination, and in particular the 
deviation from maximality, should be still regarded as quite unstable [18, 19]. 
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We have determined our best fit values minimising the quantity 


X 


2 


N 




Pi 


X?(y)-Xj 


(38) 


where the X % = X % ± ax, ’s are the N — 5 experimental parameters that we fit (in 
next Sections we will include the matter-antimatter asymmetry so that N rises to 6) 
and X\ h (y) are their predicted values depending on the theoretical parameters of the 
model y = (M[ x , M' 13 , M' 2 , M' 3 , 0; f, J) (f = 0, ±4vr/5, ±2vr/5, J=CASE A, CASE B 
are discrete parameters). Since the Dirac phase 5 can have any value within 3a, and 
has certainly not a Gaussian distribution, we decided not to include it in the fit not to 
risk to prematurely exclude potential solutions but we will comment on how this would 
change including 5. 

For the same reason for the atmospheric mixing angle we use very conservatively 
023 = 45.9° ± 3.5° both for NO and 10 (we have basically taken as a central value 
the average between the two minima and as error their halved separation). We could 
also have performed the fits distinguishing two different ranges for 0 2 3 , one in the lower 
octant for NO and one in the higher octant for 10, but these are still too weak hints. 
Our choice cuts too low or too high values for 0 2 3 , as established by all experiments 
even singularly taken, but it still treats the maximal mixing value as perfectly allowed 
and does not favour any of the two octants on the other, since, as already discussed, the 
hint coming from current experimental data is still too weak and unstable. In this way 
any emerging possible preference either for lower or for higher octant, as for maximal 
mixing, can be uniquely ascribed to the model itself and is not hidden by still unstable 
measurements. 

Within the approximation that 0 2 3 distribution is also Gaussianly distributed, the 
defined quantity y 2 has a truly ^-distribution, being the sum of squares of independent 
Gaussian variables and the p i ’ > s are the associated pulls. 

Notice that the number of (continuous) theoretical parameters is equal to the num¬ 
ber of experimental parameters (five) and, therefore the number of degrees of freedom 
vanishes. This implies that the minimisation of y 2 can only provide best fit values (and 
prediction on the low energy phases) but cannot be regarded as a goodness of fitness 
of the model, since even vanishing y 2 values could be potentially obtained with a fine 
tuned choice of the theoretical parameters y independently of the values of the experi¬ 
mental values. However, a preliminary indication of the goodness of the fitness is given 
by a comparison of the predicted value of the Dirac phase with the current best fit value 
Eq. (37), though, as we said, we certainly need a more precise measurement of 5 to draw 
firmer conclusions. 

In Table 1 and Table 2 we summarise the results of our analysis. We indicate the 
Majorana phases both in the convention diag(e* " 2 ", e* 1) and in the (PDG) convention 



CASE 

A 

£ 

+47t/5 

0 

-47t/5 

y 2 . 

A mm 

4.40 

22.8 

3.63 

mi/meV 

0.20 

0.021 

0.022 

m 2 / meV (p A m?J 

8.66 (+0.002) 

8.62 (-0.38) 

8.69 (+0.25) 

m 3 /meV (p A m?J 

49.9 (+0.76) 

48.9 (-1.34) 

49.8 (+0.47) 

#12/° ijPe 12) 

33.5 (+0.04) 

35.3 (2.37) 

33.0 (-0.66) 

$13/° (Pe 13) 

8.37 (-0.64) 

8.08 (-2.05) 

8.42 (-0.37) 

#23/° ( Pe 23 ) 

39.2 (-1.85) 

34.3 (-3.3) 

39.9 (-1.66) 

s/° 

251 

180 

109 

N° 

225 

338 

173 

N° 

82 

175 

279 

«2l/° 

217 

197 

106 

«3l/° 

135 

22 

187 

M[J GeV 

6.5 x 10 5 

6.6 x 10 5 

2.2 x 10 6 

M'J GeV 

5.0 x 10 9 

5.1 x 10 9 

5.0 x 10 9 

Mini GeV 

7.2 x 10 15 

8.6 x 10 16 

6.8 x 10 16 

M[J GeV 

3.2 x 10 10 

3.0 x 10 11 

3.4 x 10 11 

M'n/Mw 

6.4 

59.7 

69 

(j)/n 

1.27 

1.80 

1.64 

Mi/GeV 

5.1 x 10 5 

4.2 x 10 5 

4.9 x 10 5 

M 2 /GeV 

5.0 x 10 9 

5.1 x 10 9 

5.0 x 10 9 

M 3 /GeV 

7.2 x 10 15 

8.6 x 10 16 

6.8 x 10 16 


Table 1: Results for the case A (NO and no leptogenesis). 


' a 21 • a 31 

diag(l, e l , e* The table refers to NO, since we found that the model cannot repro¬ 
duce the neutrino mixing parameters for 10. In particular it badly fails in reproducing 
the neutrino mass spectrum. 

Table 1 refers to CASE A (cf. eq. (18)) while Table 2 to CASE B (cf. eq. (19)). For 
both cases we show the results for £ = ±47 t/ 5 and for £ = 0 since for £ = ±27r/5 we 
could not find any fit with a value of y 2 < 100 so that they can be basically considered 
ruled out. 

Let us list the main results postponing some comments to the conclusions. 

• We do not find any solution for 10, the model is unable to reproduce 10 neutrino 
masses. 

• For £ = ±27 t/ 5 there are no solutions with y 2 lin < 100. 

• Except for the case £ = 0, the best fit values are obtained for the CASE A. 
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CASE 

B 

£ 

+47t/5 

0 

-47t/5 

y 2 . 

A mm 

7.84 

9.32 

5.81 

mi/meV 

0.008 

3.0 

0.004 

m 2 /meV (p Aro? J 

8.58 (-0.81) 

9.25 (+0.89) 

8.72 (+0.62) 

m 3 /meV (p Am aJ 

49.6 (+0.05) 

48.8 (-1.85) 

50.3 (+1.53) 

^ 12 /° (P(9i 2 ) 

33.2 (-0.31) 

32.3 (-1.49) 

33.0 (-0.57) 

0 13/° (P0 13 ) 

8.93 (+2.10) 

8.84 (+1.66) 

8.51 (+0.06) 

^ 23 /° (P6» 23 ) 

40.0 (-1.63) 

44.55 (-0.33) 

39.9 (-1.66) 

<V° 

253 

99 

108 

N° 

335 

66 

170 

N° 

84 

325 

278 

Ct2l/° 

108 

259 

109 

«3l/° 

25 

294 

190 

M^/GeV 

7.2 x 10 5 

2.9 x 10 5 

6.6 x 10 6 

M' 22 / GeV 

4.3 x 10 10 

4.2 x 10 10 

4.4 x 10 10 

Mi/GeV 

2.2 x 10 16 

9.7 x 10 13 

4.3 x 10 16 

M[ 3 /GeV 

6.8 x 10 1U 

7.5 x 10 9 

5.1 x 10 11 

0/tt 

1.62 

1.23 

1.62 

Mi/GeV 

5.1 x 10 5 

2.9 x 10 5 

4.9 x 10 5 

M 2 /GeV 

4.3 x 10 10 

4.2 x 10 10 

4.4 x IO 10 

Ms/GeV 

2.2 x 10 16 

9.7 x 10 13 

4.3 x 10 16 


Table 2: Results for the CASE B (NO and no leptogenesis). 
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• The best fit solutions that we obtained for NO, except for the case B with £ = 0, 
seem to point to 6 2 3 ~ 40°, certainly in the first octant and even below the current 
best fit value eq. (36) hinted by current global analyses. 

• Taking into account the RH neutrino mixing parameter M' Vi improves all fits and 

one finds values of xLm( M i 3 ^ °) - = 0) - 2. 

• Values of the Dirac phase close to the best fit value from global analyses (cf. 
eq. (37)) are attained for £ = +47 t/ 5 (both in CASE A and in CASE B). If further 
experimental data should further support these values (~ —90°), this would be 
the only surviving option. 


4 Leptogenesis without flavour coupling 


Let us now consider the calculation of the matter-antimatter asymmetry with leptogene¬ 
sis. The strongly hierarchical RH neutrino mass spectrum in the A to Z model discussed 
in the previous section, with Mi <C 10 9 GeV and 10 12 GeV M 2 3> 10 9 GeV, neces¬ 
sarily points to a ^-dominated leptogenesis scenario [9] since both the lighest and the 
heaviest RH neutrino decays produce a negligible asymmetry compared to the observed 
one. In this case the B — L asymmetry, in a portion of co-moving volume containing 
one RH neutrino in ultra-relativistic thermal equilibrium, can be calculated as [9, 10] 


/V le P’ f ~ 

ly B-L ~ 


K‘ 


2e 


2 ^ 


£ 2r 0 - LK (-^2r„- L ) + £ 2e 


K. 

k 2t ± ** K( ^ 


£2 n 



+ £ 2t k(K 2t ) e s 




(39) 


where K 2t x = K 2e + K 2fl and £ 2t ± = £ 2e + £ 2fl . This expression for the asymmetry 
neglects the flavour coupling effects studied in [15]. In the present model only the last 
term will survive, as we shall justify shortly, and we may drastically approximate the 
above expression to a much simpler one, 

A T 1 b V -l — £ 2t k(K 2t ) e^ KlT , (40) 


which says that the only relevant asymmetry is that one produced at the N 2 scale 
in the tauon flavour, simply given by £ 2t k(K 2t ) with our normalisation, followed by 
exponential washout e~~^ KlT the Ni scale. 2 

2 Of course in models with Mi > 10 9 GeV or realising a crossing level solution [20] with M\ and 
close enough to have resonant leptogenesis [21], this contribution must also be considered and the 
lightest RH neutrino washout is then not a crucial problem as in the ^-dominated scenario. Examples 
of realistic S'O(IO) models realising this case were discussed in [22]. 
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The flavoured decay parameters K ia are defined as 


K = 

1 i.rv — 


r + r 

1 ia t 1 jd 

H(T = Mi 


\m Dai \- 
Miin * 


(41) 


where we introduced the neutrino Dirac mass matrix mo = v Y}' R . The Tj a ’s and the 
fj Q ’s can be regarded as the zero temperature limit of the flavoured decay rates into 
a leptons, T{Ni —> (ft l a ), and anti-leptons, T(iVj —> (f>l a ), in a three-flavoured regime, 
where lepton quantum states can be treated as an incoherent admixture of the three 
flavour components. The efficiency factors at the production can be calculated using 
[23] 


a) 


Zb{K2o) K2a 


K 2a z b( k 2a) 


z B (K 2a ) ~2 + 4Ce- 


2.5 

K 2 a 


(42) 


This expression is valid for an initial thermal abundance but, as we will see in a moment, 
since the production will prove to occur in the strong wash-out regime, there is actually 
independence of the initial RH neutrino abundance. Moreover in this case the efficiency 
factor is well approximated by K,(K 2a ) — 0.5/ 

The flavoured CP asymmetries, defined as e 2a = — (r 2a — r 2a )/(r 2 Q, + T 2a ), can be 
calculated in general as [24] 


£2 a 


3 M 2 m atm 

16 7T V 2 



(43) 


where we defined [25], 


Im 


q-OL —— 

X 2 j ~ 


m D) ia 0 m D ) aj {m ] D m D ) 


Im 


M 2 Mi m 2 m atm 


J ’OL 

2 j 
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,(44) 


with m 2 = (m' D mo) 22 /M 2 , and 


^(x) = -x 


(1 + x ) In 


1 + x 


X 


X 


X 


(45) 


Terms oc X 2l £(Mf / Mf), are strongly suppressed in a way that in the N 2 - 

dominated scenario the flavoured CP asymmetries can be approximated by 




3 M 2 m atm 
167T V 2 


2.3 


(46) 


This is because the two terms (for j = 1) from the interference with the lightest RH 
neutrinos are oc Mi/M 2 , while the second term from the interference with the heaviest 
RH neutrinos is oc M 2 /M 3 and, therefore, they are suppressed compared to Z/ 3 . 
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Let us now calculate I 23 . First of all it is immediate to see that since Y e3 = 0 then 
necessarily e 2e = 0. Even though the muon asymmetry does not exactly vanishes as 
the electronic one, however, it is suppressed as e 2 ^ oc a b 2 M[ 3 /M 3 ~ 10 -1 ' and can 
be neglected. This implies that, at least at lower order, the observed asymmetry can 
only be produced in the tauon flavour, as in Eq.(40). The expression for the final 
asymmetry, expressed in terms of the baryon-to-photon number ratio, then becomes 
extremely simple, 


Vb 


— ®sph 


nLl 

A free 
7 


^ 0.01 e 2t k(K 2t ) e T K,r , 


(47) 


where a sp h = 28/79 is the fraction of B — L asymmetry that is converted into a baryon 
asymmetry by sphaleron processes in equilibrium, k{K 2t ) is the efficiency factor for the 
tauon asymmetry at the end of the ^-production and N* ec is the number of photons 
at recombination in the given portion of co-moving volume. 3 

We have first of all to calculate K\ T and check that it is possible to have K ir < 1 . 
From the general expression eq. (41) one has 


_ \m D Ti\- 

Mi m* 


v 2 a 2 


m* Mi 


ii-7i 


(48) 


showing that the quantity 7 , dehned in Eq. (16) and originating from the mixing pa¬ 
rameter M[ 3 , plays a crucial role. Indeed since v 2 a 2 /(m*Mi) ~ 20 S> 1, the possibility 
to have K\ T < 1 necessarily relies on having Arg[y] ~ 0 implying q i rs_/ -(£ + 3tt/5). 

Let us now calculate e 2r from the eq. (46). Considering that m 2 = 21 v 2 b 2 /M 2 and 
that ZJ 3 ~ 4 b 2 c 2 sin(2 0 + £ + 6 7 t/5 ) v 4 /(M 2 M 3 m 2 m atm ) , one hnds 

c 2 M 

£ 2 t - 7 ^- -Tj- sin(20 + i + 6tt/5) . (49) 

Zo IT 1V1 3 

When the condition for K\ r < 1 on 0 is imposed one has 


£2 T 




c 2 M 2 
28 it M 3 


sin , 


(50) 


showing that out of the five possible values of £, only ^ = —27r/5, — 47 t /5 can lead 
to the correct sign of the asymmetry. An order-of-magnitude estimation gives then 
|e 2T | ~ 10- 7 -10" 6 . 

Finally from the eq. (41) we can calculate 


\m DT2 \ 2 4 b 2 v 2 
M 2 m * (9) M 2 m* 


(51) 


3 This expression is valid independently of the normalisation of the abundances Nx, since the nor¬ 
malization factor would cancel out in the ratio N i B _ L /NZ ) ec . 
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that, plugged into the eq. (47) gives t]b ~ 10 _ 11 -10 -9 , showing that potentially the 
observed value of the asymmetry could be reproduced. However when one tries to fit 
simultaneously the asymmetry and the mixing parameters one finds that the condition 
7 ~ 1, necessary to have K\ T < 1, is incompatible with the possibility to reproduce the 
correct values of the mixing parameters so that the asymmetry produced at the N 2 scale 
is afterwards completely washed out at the N\ scale. 


5 Leptogenesis with flavour coupling 


As we have just seen the asymmetry is mainly produced by the next-to-lightest RH 
neutrinos in the tauon flavour but this asymmetry is fully washed-out by the lightest 
RH neutrinos since the condition K lr < 1 is not compatible with the measured values 
of the mixing parameters. 

However, one has also to consider that part of the asymmetry in the tauon flavour is 
transferred to the electron and muon flavours by flavour coupling effects due primarily 
to the fact that V 2 -decays produce in addition to an asymmetry in the tauon lepton 
doublets also an (hyper charge) asymmetry in the Higgs bosons. This Higgs asymmetry 
unavoidably induces, through the inverse decays, also an asymmetry in the lepton dou¬ 
blets that at the production are a coherent admixture of electron and muon components. 
Therefore, in this case, inverse decays actually produce an asymmetry instead of wash it 
out as in a traditional picture. A somehow smaller effect is also due to the asymmetries 
stored into quarks and into right handed charged leptons. 

The account of flavour coupling effects modifies the usually considered expression 
for the asymmetry eq. (47) resulting into [15] 

Vb ^ > (52) 

ct=e,/x,r 


where in addition to the tauon contribution eq. (47) one also has an electron contribution 


V ( b - -0.01 e 2 r k(K 2t ) 


I<2e 


K 2e + 

and a muon contribution given by 


C {2) e~^ Kle 


(53) 


r] { B - —0.01 £ 2t k(K 2t ) 


K‘ 


2/x 


K 2e + K, 


C 


( 2 ) 


2[i 


jdy W3) 
Klr ^ 


3tt is 
0 8 


(54) 


It should be noticed how the source of the electron and muon asymmetries is in any case 
the tauon asymmetry, but part of this induces a muon and an electron asymmetry thanks 
to flavour coupling. The flavour coupling coefficients are given by C^ 2 } = 104/589 and 
cjir = 142/537. Notice that we are neglecting additional correcting terms containing 
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products of the flavour coupling coefficients and we are also neglecting terms oc e 2r ± 
since this is always too small to generate sizeable contributions. Also notice that from 
Eqs. (15) and (41) one can see immediately that K j e = 0. 

Notice moreover that the second term in the muon contribution comes from flavour 
coupling at the lightest RH neutrino wash-out that works in the same way as at the 
production: the lightest RH neutrino wash-out processes acting on muon lepton dou¬ 
blets, in the presence of a non-vanishing Higgs asymmetry, induce a muon asymmetry. 
However, this term is basically much smaller than the first term that gives the dominant 
contribution to the total asymmetry, indeed it dominates on the electronic term as well 
since K 2e -C A" 2/t , and is in the end responsible for the two solutions that we found. 

We have indeed performed again a numerical fit for all the different cases and this 
time we have found that, both for the case A and for the case B with £ = + 47 t/ 5 , there 
is indeed a solution, shown in Table 3, with an acceptably small y 2 value. We found no 
other solution with a y 2 lower than 100 for all the other cases (£ = 0, ±27 t/5, — 47 t/ 5 ) 
since these tend to give either a too small 9 32 or a too small $13 or both. In both cases 
one can see that K \ /t <1. In this way the flavour coupling induced asymmetry at the 
production in the muon flavour can survive the lightest RH neutrino wash out giving the 
dominant contribution to the final asymmetry. This is a nice example of the “flavour 
swap” scenario envisaged in [15]. 

It should be noticed that both solutions predict the atmospheric angle well in the 
second octant, a feature that will be relatively soon tested by new data from neutrino 
oscillation experiments. In fact the two solutions give quite similar predictions on 5 
and $ 23 - hi the Figure we have plotted, in the plane 5 vs. 9 23 , points corresponding to 
solutions with y 2 < 10 . This gives an idea of the allowed region in this plane. It can be 
seen how for the CASE A the minimum values correspond to (S,9 23 ) ~ (10°, 51°) while 
for the CASE B the whole region is slightly reduced (indeed Xmin — 5 for the CASE A 
and Xmin — 6 for the CASE B) and shifted to higher values both of 5 and 9 23 and the 
minimum values are (S,9 23 ) ~ (14°, 52°). This shows that the CASE A is slightly more 
favoured compared to CASE B. 

It should also be noticed that since now the number of degrees of freedom v — 1, the 
Xmin can be regarded as an indication of the goodness of fit (g.o.f.), having in mind the 
previous discussion on the measurement of 9 23 . Of course we also had the possibility 
to choose the value of the discrete parameter £ and between CASE A and CASE B 
and this somehow made things a bit easier, but it is still intriguing that, despite its 
reduced number of parameters, the model can also account for the matter-antimatter of 
the Universe. It should also be said that if the current experimental information on 6 
is taken into account (cf. eq.(37)), one probably would have an additional contribution 
Ax 2 = p 2 > 3 so that the allowed regions shown in the Figure are marginally compatible 
with current data on <5, especially in CASE B while in CASE A the portion around the 
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CASE 

A B 

£ 

+47t/5 

Xmin 

5.15 

6.1 

M( 1 /10 6 GeV 

1.30 

1.33 

M' 2 /10 10 GeV 

0.48 

4.35 

M' 3 /10 12 GeV 

2.16 

1.31 

M( 3 /10 10 GeV 

1.81 

0.61 


3.75 

0.141 

0 /tt 

0.795 

0.788 

M^IO 7 GeV 

15 

2.7 

M 2 /10 10 GeV 

0.483 

4.35 

M 3 /10 12 GeV 

2.16 

1.31 

N 

203 

38 

mi/meV 

2.3 

2.3 

m 2 /meV (p A m?J 

8.93 (-0.22) 

8.94 (-0.25) 

m 3 /meV (p Am aJ 

49.7 (+0.17) 

49.7 (+0.21) 

St m i/ meV 

61 

61 

m ee /meV 

1.95 

1.95 

^12/° (p<? 12 ) 

33.0 (-0.58) 

33.0 (-0.66) 

0is/° (P0is) 

8.40 (-0.47) 

8.40 (-0.49) 

023/° (pe» 23 ) 

53.3 (+2.1) 

54.0 (+2.3) 

<V° 

20.8 

23.5 

N° 

118 

115 

N° 

281 

278 

Oi2\/° 

163 

162 

«3l/° 

242 

245 

P/° 

279 

279 

o!° 

220 

221 

W io ~ 10 (zO 

6.101 (+0.01) 

6.101 (+0.01) 

^2r 

-8.1 x 10" 6 

-1.3 x 10" 5 

/+„ 

0.11 

0.58 

K i T 

4341 

800 

I<2r 

7.3 

7.3 

k 2ii 

29.2 

29.2 

I<2e 

1.8 

1.8 

|(m !2 ) 11 |/meV 

6.6 x 10“ 6 

3.7 x 10" 5 

|(m“ 1 )3 3 |/meV -1 

0.22 

1.2 

|(m !2 ) 12 |/meV 

2.6 x 10" 5 

1.5 x 10~ 4 


Table 3: Solutions found for flavour coupled leptogenesis with xLin < 100 (neutrino masses are NO). 
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Figure 1: Scatter plots of points in the plane S vs. 6*23 with y 2 < 10 (leptogenesis included) for the 
CASE A (left) and for the CASE B (right) for NO and £ = +47 t/5 . The red points correspond to the 
best fits solutions in Table 3. These predictions are subject to the theoretical uncertainties discussed 
in the text. 


minimum values for (<5, $ 23 ) is still allowed at ~ 2 a. 1 However, one should also take into 
account that we are currently neglecting corrections from the charge lepton mass matrix 
that has been approximated to be diagonal and also uncertainties on the values of the 
up quark masses at high scale (see Footnote 1). This might help in shifting the allowed 
regions shown in the Figure toward more favoured experimental values. In any case the 
allowed regions shown in the Figure in the plane 5 vs. 0 23 , in combination with NO, are 
quite a strong prediction of the model that will be certainly tested during next years. 
Much more difficultly testable predictions, certainly not in a close future, are the very 
small value of the neutrinoless double beta decay effective neutrino mass ( m ee ~ 2 meV) 
and the small increase (~ 3 meV) of the sum of the neutrino masses from its hierarchial 
value {J2i m i\ mi = 0 — 58meV). 

A significant feature of the best fit spectrum for both cases A and B in Table 3 is 
the compressed spectrum of right-handed neutrino masses as compared to the original 
estimates in Eq.(4). This stems from the requirement that successful N 2 leptogenesis is 
proportional to the ratio M 2 /M 3 in Eq.(50) which must necessarily be larger than naively 
expected in Eq.(4). M 2 cannot change much, since according to sequential dominance 
and the charm quark mass relation in Eq.(26), it sets the solar neutrino mass scale. It 
follows that the only possibility is for M 3 to decrease which requires the parameter m c 
to increase (see Eq.(26)) and hence the third seesaw matrix in Eq.(23), proportional to 

4 However, note that in Eq.(38) we would then have N = 7 so that the number of degrees of freedom 
would be 7-5=2 so that one could say that including S would actually improve the fit, decreasing 
y 2 /d.o.f, for points with Ay 2 < 5. This point is right now quite indicative since the distribution of 5 
is highly non Gaussian. 
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diag(0, 0,1), threatens to dominate the neutrino mass matrix and pull the atmospheric 
mixing away from its maximal value. This threat is averted by noting that, for large 
7 , the first matrix in Eq.(23) has a dominant (3,3) element of order y 2 which may 
partially cancel the contribution from the third matrix leaving a resulting (3,3) element 
of order 7 , of the same order as the (2,3) and (3,2) elements of the first matrix. This 
requires a fine tuning of one part in 7 (tabulated in Table 3). In order to maintain the 
correct atmospheric neutrino mass, this increase in 7 must be compensated by reducing 
the parameter m a in proportion to 7 which, from Eq.(26), requires My to increase in 
proportion to 7 . In the next section we obtain further analytic insight into our results 
from different perspectives. 

6 Analytical insight from general S'O(10)-inspired lep- 
togenesis 

It is interesting to get some understanding of the numerical results we obtained within 
the A to Z model from the more general context of S'O(10)-inspired leptogenesis [26, 
20, 11, 12]. In this case the asymmetry can be calculated analytically assuming that 
the spectrum of the neutrino Dirac mass matrix is hierarchical {rn\\ 3> m® 2 3> nr^), 
an assumption certainly holding in our case. Here we present a simple generalisation of 
the results presented in [12] providing a very good explanation of the numerical results. 

The starting point is to write the neutrino Dirac mass matrix in the bi-unitary pa- 
rameterisation (mathematically equivalent to the singular value decomposition of m^), 

mD v = Vl D m° U R , (55) 

where we defined D m D = diag(m^, m^ 2 , m® 3 ). The unitary matrices Vl and Ur trans¬ 
form respectively the LH and RH neutrino fields from the flavour basis to the Yukawa 
basis. In our case, from the eq. (15) for the neutrino Yukawa matrix and parameterising 
Vl analogously to the leptonic mixing matrix as 

Vl = R2-M,) RM 3 ) R 12 (9 l 12 ) D($ l ) , (56) 

one has tiny 9 23 ,9f 3 <C 1° so that Vl — Ri 2 (9f 2 ) with 9 32 ~ 14°. Notice that this 
angle would correspond to the Cabibbo angle 9c in the quark sector, that is therefore 
overestimated by ~ 1°. However, turning on down quark mass matrix (non-diagonal) 
terms, one can reproduce the correct value. Similar corrections are expected on the 
neutrino mixing angles from analogous charged lepton mass matrix correcting (non- 
diagonal) terms. 

Inserting the eq. (55) into the see-saw formula eq. (22) written in the flavour basis, 
one obtains the following expression for the Majorana mass matrix in the Yukawa basis 

Mr = U R D m U ] r = -D m u Vt m~ l v£ D m n , (57) 
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and for its inverse 


U R Djf U\ = —-D“d V l m v V[ D~^ . (58) 

In the above we have transformed the heavy Majorana mass matrix from the diagonal 
basis Mr = Dm in Eq.(ll) to the basis in Eq.(55) in which the Dirac neutrino mass 
matrix is diagonal. The analytical expressions for the RH neutrino masses obtained in 
[12] in the approximation Vl — I get extended for a general Vl simply replacing m u —y 
m v = Vl m u V^ [20] (the light neutrino mass matrix in the Yukawa basis) obtaining for 
the RH neutrino masses the following analytical expressions 

\{fh v ) u | 


Mi ~ 


«i) 2 


(m v In 


M 2 ~ 


« 2 ) 2 


M 3 - « 3 ) 2 IK/)33 


mi m 2 m 3 | (m 




/33 


(59) 


In deriving the first and third equalities above we have used the strong up-type quark 
mass hierarchy (equal to the eigenvalues of the neutrino Dirac mass matrix in case A), 
and for the second equality we have taken the absolute value of the determinant of the 
seesaw formula. Plugging the expressions for M 2 and M 3 into the eq. (49) for e 2r and in 
turn this into the eq. (54) for the dominant muonic contribution to Vb and taking into 
account that K 2fl /(K 2fl + K 2e ) = 16/17, one finds 

0.04 c 2 M 2 

Id 


Vb - 


- -k(K 2t )Cv e 

119 7T M 3 v ; T±T 


( 2 ) 0 -^-K 


8 sin(2</> + £ + 67 t/5). 
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From the expressions for M 2 and M 3 in eq. (59) and, taking into account the rela¬ 
tions (17), (20) and (21), one can also write for CASE A (CASE B) 

0.04 (9)m 2 harm |(m l ,) 11 | 
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where we used the approximations sin 9c cos 9c ~ 1 and the following analytical 

expression for Ur [20, 12] 
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that we generalised here to the case when Vl ^ I and where the expressions for the 
phases given in [12] can be also generalised in terms of m u (we do not need them 
here). From these expressions for we can now make some considerations that explain 
some of the features of the two found numerical solutions for £ = +47 t/ 5 (see Table 3) 
providing a useful analytical insight. 
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• Inserting the numerical values for M 2 /M 3 , Jl lM and 0 given in Table 3 and 4 into 
the eq. (60), the observed value of r/ B is indeed reproduced. 

• For one of the possible choices of values for £, the effective lcptogenesis phase is 
maximal for cj) = 37t/20 — £/2 + n n. For the case £ = 47t/ 5 the phase is maximal 
for cf) = 37t/4 + rm explaining quite well the best fit values found for 0 (see Table 
3). 

• From the eqs. (59) we can derive an analytical expression for M 2 /M 3 , the other 
crucial parameter determining the value of the asymmetry (cf. eq. (60)), hireling 

^2 = (mg 2 ) 2 \(m u ) n \ (61) 

M 3 (m^ 3 ) 2 m i m 2 m 3 Km” 1 )^! 2 ’ 

We have verihed that indeed, inserting into this equation the measured values for 
the mixing angles and for the solar and atmospheric neutrino mass scales, one 
obtains the correct value of M 2 /M 3 thanks to a reduction of |(m“ 1 ) 33 | to val¬ 
ues ~ O.leV -1 from phase cancellations (without cancellations one would have 
|(fh“ 1 ) Tr | ~ 100eV _1 ) and that this is necessarily accompanied by the reduction 
of \(m u )u\ to values ~ 10 -6 eV (without cancellations one would have |(m.^)n| ~ 
10 nreV) when a condition K\ tl < 24 is also imposed (it is then much more general 
than Kin ~ !)• One can see th a t this is indeed what happens from Table 3, where 
we also show the best fit values for |(m ! 7 1 )33|, |(uA)n| and |(m iy )i 2 |. It is interest¬ 
ing that this results into a stable value of the next-to-lightest RH neutrino mass 
(M 2 ~ 5 x 10 9 GeV for the CASE A and M 2 — 2 x 10 10 GeV for the case B) the 
same we obtained in Table 2 without imposing lcptogenesis, while the heaviest RH 
neutrino mass reduces to values M 3 ~ 2 x 10 12 GeV. Correspondingly the lightest 
RH neutrino mass, though it does not play a direct role since the asymmetry is 
V 2 -dominated, is necessarily forced to grow to values M, rsj 10 7 ^ 8 GeV. This is 
something we have already discussed at the end of Section 5 starting within the 
model parameterization and that we have now seen also from a bottom up per¬ 
spective. Moreover the value of the lightest neutrino mass and of the neutrinoless 
double beta decay effective neutrino mass mi ~ m ee ~ 2 meV, in agreement with 
what we obtained fully numerically in Table 3. 

Therefore, the eqs. (59) do explain quite well the obtained results and they also 
seem to point to a more general result: given the measured values of the low en¬ 
ergy neutrino parameters, the reduction of M 3 generally (the condition Ah /t < 24 
is quite general) implies an uplift of Mi while M 2 remains stable at a scale 
~ 10 10 GeV. In other words the current measurements are such that the con¬ 
dition for the realisation of the M 2 — M 3 crossing level (reduction of | (m.^ 1 ) 33 1) 
necessarily implies also the occurrence of the Mi — M 2 crossing level solution, that 
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would lead to a quasi-degenerate RH neutrino mass spectrum in the close vicinity 
of the crossing level. We have checked this statement also verifying that if one let 
#13 to be free, then this effect occurs only for # 13 ~ 7° 4- 15°. In other words, with 
current values of the neutrino oscillation parameters, a deviation from a strong 
hierarchy of the RH neutrino masses proportional to the squares of the up-quark 
masses, seems to point, for K\ tl < 24, toward a RH neutrino spectrum where all 
three RH neutrino masses tend to get closer. This analytical insight seems to help 
understanding also recent numerical results where quark-lepton parameters have 
been fitted within S'O(10)-models either excluding [27] or including [28] lcptogen- 
esis and either hierarchical or compact RH neutrino mass spectra have been found 
but never crossing level solutions with only two close RH neutrino masses. 

• The best fit values of the flavoured CP asymmetry £ 2 T in Table 3, the source of 
the asymmetry, is well above the upper bound [ 8 , 29, 25] 


3 Mi m atm IK\ a _ 6 / Mi \ /K la 

16 vr v 2 V Ki ~ Vl0 10 GeVy V ’ 


(65) 


holding for the term in the lightest RH neutrino CP flavoured asymmetries anal¬ 
ogous to the first term in eq. (43). Indeed one has |2J 3 | ~ 17 3> 1 for the best fits 
of Table 3. This is possible if the absolute values of the entries of the orthogonal 
matrix are much larger than unity. This can be understood considering that the 
orthogonal matrix [30] in 5O(10)-inspired models is given by (U = Upmns ) 
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where = diag(e 'ye , e ^a 3 ), an analytical expression that generalises 
that one given in [12] when Vj j ^ I. One can see that for the best fit values of rhn 
and mjg 1 one has that all the absolute values of the f! entries in the first and third 
column are much higher than unity. This is confirmed by the numerical results we 
find for Q, corresponding to the best fits in Table 3 for CASE A 


^(CASEA) _ 


-4.40016 - 15.9889 i 0.0930875 - 0.894045 i 
-15.9446 + 3.40333 i -1.15394 + 0.0537137i 
-3.69174 + 4.35811 i 0.709793 + 0.204576 i 


— 16.0396 + 4.38107 i 
3.40494 + 15.9553 i 
4.37787 + 3.64191 i 


.(67) 
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and for CASE B 


/ -1.77835 -6.85986 * 0.108413 - 0.897431 * -6.97828 + 1.73423 * \ 

q(caseb) I -6.87598 + 1.34103 i -1.15331 + 0.0386159 * 1.34278 + 6.90018* 

\ -1.64314 + 1.81259* 0.710523 + 0.199612* 1.85785 + 1.52677* 

These specific expressions for 11 should give a clear idea of the involved fine 
tuned cancellations and explains why one can have |XJ 3 | ;+ 1 (implying e 2r 
(3M 2 m atm /(16 Tiv 2 )) and more generally why the flavoured CP asymmetries can 
be enhanced in the vicinity of crossing level solutions [20] though RH neutrino 
masses are still hierarchical and leptogenesis is far from being resonant. However, 
our novel solution relying on flavour coupling, represents in this respect quite 
a large improvement compared to the commonly considered solutions based on 
M ! > 10 9 GeV, since this requires a further uplift of the lightest RH neutrino 
mass of at least two orders of magnitude and, therefore, even higher fine tuned 
cancellations. 

• As we have seen the third important ingredient for the existence of the solution is 
to have K i jx < 1 . This is crucial in the case A, while in the case B the electron 
contribution, sixteen times smaller than the muonic one, would be still sufficient 
(remember that K\ e = 0). The expression for K\ tl eq. (62) shows that this condi¬ 
tion is realised for |(m I/ )i 2 | 2 < m* | (rri^)n | that is indeed verified for the best fits 
in Table 3. 

• We can compare our results to those presented in [11] for Vl — Vckm where flavour 
coupling was neglected. The values of the Majorana phases in Table 3, match with 
those in [11] though in a marginal region (for an easier comparison in Table 3 we 
also give the the Majorana phases in the convention diag(e ip , 1, e la ) as in [20, 11]). 
The same it is true for the values of 6 13 , # 23 , 6 . In [11] it can be also seen how the 
ratio M 2 /M 3 can get reduced around m\ ~ 2.5 meV. Moreover in [11] there is a 
lower bound on M 2 > 5 x 10 10 GeV that corresponds to a 2 = m^ 2 /m c h arm > 3, 
while we find that successful leptogenesis is obtained for M 2 ~ 5 x 10 9 GeV and 
we have a 2 = 1 in the CASE A. This is because the CP asymmetry, thanks to 
the condition K\ fl that makes possible much lower values of M 2 /M 3 , is greatly 
enhanced compared to the upper bound eq. (65). Moreover in [11] there are no 
points with M\ uplift. This is explained since there the solutions correspond to 
K\ r < 1 while here with flavour coupling the solutions correspond to A' 1/t < 1 
and so there is an intrinsic difference. This shows that the account of flavour 
coupling indeed opens new solutions enlarging the allowed regions in the space of 
parameters, in particular making possible a reduction in the scale of leptogenesis 
set by M 2 . 


.(68) 
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7 Conclusions 


The A to Z model can not only provide a satisfactory fit to all parameters in the 
leptonic mixing matrix but can also reproduce the correct value of the matter-antimatter 
asymmetry with A 2 -dominated leptogenesis. In this respect it is crucial to account for 
flavour coupling effects due to the redistribution of the asymmetry in particles that 
do not participate directly to the generation of the asymmetry, in primis the Higgs 
asymmetry. In particular a “flavour swap” scenario is realised whereby the asymmetry 
generated in the tauon flavour emerges as a surviving asymmetry dominantly in the 
muon flavour. The solution works even in the simplest case where the neutrino Dirac 
mass matrix is equal to the up quark mass matrix. 

Neutrino masses are predicted to be NO, with an atmospheric neutrino mixing angle 
well into the second octant and the Dirac phase 5 ~ 20°, a set of predictions that will be 
tested in the next years in neutrino oscillation experiments. We expect these values to 
be slightly corrected by charged lepton mass matrix corrections and different theoretical 
uncertainties (for example in the values of the up quark masses at the high scale). In 
particular we note that charged lepton mixing corrections, although small in the A to 
Z model due to the (1,2) entry of charged lepton and down quark mass matrices being 
zero, may yield atmospheric mixing corrections of order one degree. 

In conclusion, the novel solution that we presented involving indispensible flavour 
coupling, opens new possibilities for successful leptogenesis within realistic 50(10)- 
inspired models. Although there is fine tuning given by the parameter 7 , the level of 
fine tuning is mild, certainly much smaller than in traditional quasi-degenerate solu¬ 
tions. The reason for the fine tuning is that the spectrum of RH neutrinos, must be 
compressed as compared to the “natural” 5O(10)-inspired spectrum proportional to the 
squares of up-type quark masses. However, whereas models with M\ > 10 9 GeV, com¬ 
monly considered in the literature, involve a very compact mass spectrum, approaching 
degeneracy, the compressed mass spectrum for the W-dominated case considered here 
remains very hierarchical. 
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